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ABSTRACT 

The origin of grain size distribution in the interstellar medium is one of the most fundamental 
problems in the interstellar physics. In the Milky Way, smaller grains are more abundant 
in number, but their origins are not necessarily specified and quantified. One of the most 
efficient drivers of small grain production is interstellar turbulence, in which dust grains can 
acquire relative velocities large enough to be shattered. Applying the framework of shattering 
developed in previous papers, we show that small (a < 0.01 |^m) grains reach the abundance 
level observed in the Milky Way in ~ 10 s yr (i.e. within the grain lifetime) by shattering in 
warm neutral medium. We also show that if part of grains experience additional shattering in 
warm ionized medium, carbonaceous grains with a ~ 0.01 |xm are redistributed into smaller 
sizes. This could explain the relative enhancement of very small carbonaceous grains with 
a ~ 3-100 A. Our theory also explains the ubiquitous association between large grains and 
very small grains naturally. Some tests for our theory are proposed in terms of the metallicity 
dependence. 
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1 INTRODUCTION 

Formation and evolution of dust grains is one of the fundamen- 
tal problems in the interstellar physics. In particular, the grain 
size distribution largely affects the opacity (or extinct i on) o f inter- 
stellar medium (ISM). iMathis, Rumpl. & Nordsieckl dl977l . here- 
after MRN) reproduced the extinction curve of the Milky Way 
by a mixture of silicate and graphite with a size distribution of 

■ n(a) oc a -3,5 , where n(a) da is the numbe r density of grains with 
radii b etwe en a and a + da. M ore recently, IWeingartner & Draind 

\ d200lh and iLi & Draine] feOOlh have derived more detailed grain 
size distribution, finding that an enhancement of small carbona- 
ceous grains, especially polycyclic aromatic hydrocarbons (PAHs), 
is required to reproduce the mid-infrared (MIR) spectrum in the 
Milky Way. 

The grain size distribution reflects the physical proce sses in 
grain formation and evolution dO'Donneir & Mathis 1997). Dust 
grains are not only produced and ejec ted by supernovae (SNe) (e.g. 
iKozasa. Hasegawa. & Nomotdll989|) and asympt otic giant branch 
(AGB1 stars (e.g. [Gail. Keller. & Sedfmayrl 1 1984) but also pro- 
cessed in the ISM. Hirashi ta & Yanl d2009T hereafter HY09) show 
that shattering occurs efficiently in diffuse ISM, where grains are 
accelerated by magnetohydrodynamic (MHD) turbulence. In par- 
ticular, turbulence in warm ionized medium (WIM) can drive shat- 
tering o f dust grains in < 1 Myr, producing small grains effi- 
ciently I Hirashit a et al J[201oh . Therefore, shattering should play a 
significant role in determining the grain size distribution. Shattering 
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also occurs efficiently in supernova sho cks dBorkowski & Dwekl 
ll995l : |jones. Tielens. & Hollenbacldl 199(1 . 

There are some indications that the gr ains formed and sup- 
plied to the ISM are biased to large sizes. As lNozawa et al. 1 1200% 
show, the grains ejected from SNe II are biased to large sizes be- 
cause they suffer from the destructio n in the shocked region be - 
fore being ejected into the ISM (cf. iBianchi & Schneider] 120071) . 
The grain radius is typically around a ~ 0.1 \im. The typical size 
of grains condensed in AGB stars is also suggested to be large 
(a ~ 0.1 urn) from the observations of sp ectral ener gy distributions 
Groenewegenll 19971 ; iGauger et al.lll999h . although Hofmann et al.l 
200 lh show that the grains are not single-sized. Moreover, in the 
Milky Way and other systems whose metallicity is around solar, 
the major part of the dust mass comes from the gra in growth b 
accretion of heavy elements in interstellar clouds (e.g. lDwe " 
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lDrainell2009h . Not only the grain growth by accretion but also co- 
agulation in dense clouds causes a strong depletion of small grains 
(HY09). To summarize, the grains supplied from molecular clouds 
or stars should be biased to the largest size range, a ~ 0.1 \im. 

If the major part of dust grains supplied are large, the origin of 
small grains is worth investigating seriously. Here we focus on shat- 
tering as a production source of small grains. Specifically, we ex- 
amine if large grains with a ~ 0.1 (im are efficiently redistributed 
into small (< 0.01 irm) grains by shatterin g driven by interstellar 
turbulence, which is ubiquitous in the ISM. lHirashita et alj feOlOh 
have shown that shattering in the WIM produced by starburst activ- 
ities can efficiently produce small grains from large grains supplied 
from SNe II. In this Letter, we discuss shattering in a more gen- 
eral context by showing that turbulence plays a fundamental role in 
supplying small grains in more quiescent environments such as in 
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warm neutral medium (WNM) covering a large fraction of the ISM 
in the Milk y Way. Shattering in supernova shocks should also be 
considered \ Jones et alJl996h . but the inclusion of this process into 
our framework is left for future work. Throughout this paper, grains 
are assumed to be spherical, and the words, 'large' and 'small', are 
used for grain radii of a ~ 0.1 |im (~ the largest size range in 
MRN) and a < 0.01 pin, respectively. 



HY09 consider shattering in various ISM phases. Among 
them, we first consider WNM, which occupies a significant frac- 
tion of the interstellar space. We also treat additional shattering in 
WIM. Although WIM occupies a small fraction in the ISM, the 
short time-scale of shattering in WIM (HY09) can make a signifi- 
cant imprint on the grain size distribution. 



2.3 Grain velocities 



2 METHOD 

2.1 Initial condition 

To examine interstellar shattering induced by turbulence as a mech- 
anism of producing small grains from large grains, we assume 
the initial grain size distribution to be dominated by large grains. 
This 'initial' grain size distribution effectively represents not only 
one at the formation in stellar ejecta but also one after the grain 
growth in cold clouds in the ISM. The functional form of the initial 
grain size distribution is represented by a log-normal distribution, 
n(a) = n ini (a): 



C 

ni n i(a) = { a 




exp 



ln(a/a )] 2 
2cr 2 



(otherwise), 



(i) 



where G is the normalizing constant determined by equation lO, 
ao and a are the central grain radius and the standard deviation of 
the log-normal distribution, respectively, and a m j n and a max are 
the minimum and the maximum grain sizes, respectively. We adopt 
ao = 0.1 pm according to the above argument in Introduction, 
own = 3.5 A, a m ax = 1 Vim, and a — 0.6. The selection of a is 
ba sed on the size dispe rsion of large Si and C grains as calculated 
bv lNozawa et al.1 d2007h . 

We consider two dust species, silicate and graphite, and dis- 
tribute each species according to the above size distribution. The 
normalizing constant G is determined for each species by 



4yr 3 / \ j 

—a pg r n in i(o) da, 



(2) 



where 1Z is the dust-to-hydrogen mass ratio (i.e. dust abundance 
relative to hydrogen), nu is the hydrogen number density, and p gI 
is the grain material density (3.3 g cm" 3 and 2.2 g cm" 3 for silicate 
and graphite, respectively). We adopt 1Z = 4.0 x 10~ 3 and 3.4 x 
10~ 3 for silicate and graphite, respectively, according to the typical 
Galactic dust-to-gas ratio (HY09). 



is calcu- 



2.2 Shattering 

The evolution of grain size distribution by shatterin 
lated based o n HY09, whose formulation is taken from Jones et all 
dl994lll996h . The shattering equation is calculated for silicate and 
graphite separately to avoid the complexity caused by the colli- 
sion between different species. The grain-grain collision rate is 
estimated based on the grain velocity as a function of grain size 
(Section [2.3t , and if the relative velocity is higher than the shatter- 
ing threshold (2.7 and 1.2 km s~ for silicate and graphite, respec- 
tively), the grains are redistributed into smaller grains according to 
the size distribution of shattered fragments, which is assumed to be 
power-law (o c a~ 3 3 ; note that the results are not ve ry sensitive to 
the exponent; I Jones et al. 1996; Hirash ita et alj201(il) . 



The grain veloci ty as a function of grain radius a is taken from 
lYan et al.1 d2004l) . who considered the grain acceleration by hy- 
drodrag and gyroresonance (resonance which occurs when the 
Doppler-shifted frequency of the MHD wave in the grain's guid- 
ing centre rest f r ame i s a multiple of the gyrofrequency) based on 
lYan & Lazarianl ((2003) and calculated the grain velocities achieved 
in various ISM phases. They adopted the following physical con- 
ditions: in WNM n H = 0.3 cm" 3 , T = 6000 K, n c = 0.03 
cm" 3 , Guv = 1, Va = 20 km s" 1 , L = 100 pc (T is the gas 
temperature, n c is the electron number density, Guv is the UV 
radiation field normalized to the typical Galactic value in the so- 
lar neighbourhood, Va is the Alfven speed, and L is the injection 
scale for turbulence), and in WIM n H = 0.1 cm" 3 , T = 8000 
K, n — 0.0991, and the same values as in WNM for the other 
quantities. 

Large (a > 0.1 (am) grains are mainly accelerated by gyrores- 
onance. Thus, in considering shattering of large grains, gyroreso- 
nance is the key process. Although the grain charge strongly de- 
pends on gas density, electron fraction, and UV flux, the veloci- 
ties of large grains do not strongly d epend on the grain charge as 
mentioned in Hirashi ta & Yanl ( 120091) : The acceleration by gyrores- 
onance increases with the grain charge, but the acceleration du- 
ration (the hydrodrag t ime-scale) decreases with the grain charge 
jYan & Lazarianl 120031) . Thus, as for the small grain production 
from large (a > 0.1 [im) grains, the re sults in this paper is insen- 
sitive to the grain charge. As shown in iHirashita et al.l ( |2010|) . the 
velocities of large grains do not sensitively depend on the gas den- 
sity. The grain velocity may depend on the electron fraction, since 
the damping of turbulent motion depends on it. The dependence 
on the ionization fraction is roughly bracketed by the velocities in 
WNM and WIM. Since the velocities of large grains are similar in 
these two phases, they are not affected by the electron fraction. The 
remaining uncertainty in the model comes from the magnetic field 
(or the Alfven velocity), which is discussed in Section [3~T1 



3 RESULTS 

3.1 Shattering in WNM 

The evolution of grain size distribution by shattering in WNM 
is shown in Fig. Q] As a reference, the MRN size distribution 
(n(a) oc a" 3 ' 5 with a m i n = 10~ 3 pin and a max = 0.25 ^m) is 
also shown. We observe that shattering redistributes the large grains 
into small sizes on a time-scale of < 100 Myr, wh ich is shorter tha n 
the grain lifetime in the ISM (a few xlO 8 vr: Ijones et al.lll996h . 
Moreover, the grain size distribution approaches the MRN size dis- 
tribution, which indicates that shattering is efficient enough to en- 
rich the small grains up to the abundance level in the Milky Way. 

Since the time-scale of shattering is inversely proportional to 
the gas density and the dust-to-gas ratio, we can summarize our 
result by the following estimate of a time-scale, r s h a t, on which 
shattering in turbulence supplies small grains to a level compatible 
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Figure 1. Grain size distributions in WNM. The grain size distribution is divided by njj and multiplied by a 4 to show the grain mass distribution per log a and 
per hydrogen nucleus. Two grain species, (a) silicate and (b) graphite, are shown. The solid and dot-dashed lines indicate the grain size distribution at t = 50 
and 100 Myr, respectively, and the dashed line shows the initial log-normal distribution. The dotted line shows the MRN size distribution. 



with the MRN size distribution: 

x (ao/0.1 nm)(V shat /20 km s" 1 )'" yr, (3) 

where TZq is the dust-to-hydrogen mass ratio assumed in this paper 
for the Galactic environment, and V s hat is the typical relative ve- 
locity among the grains. The dependence on ao comes from the 
grain collision time-scale with a fixed total grain mass. For the 
grain velocity, V s ha.t ~ Va can be achieved after gyroresonance. 
The exponent a is 1 if small grain production is only governed by 
the collision rate. We test several cases where we have artificially 
changed the velocities achieved by large grains, and a proves to be 
~ 2 around Khat ~ 5-30 km s _1 . As mentioned in Section l4~2l 30 
per cent of the small grains are produced by the collision between 
small and large grains; in this case, cratering occurs in large grains, 
and a becomes larger than 1 because of the velocity dependence of 
the cratered volume. 

Shattering of large grains in cold neutral medium (CNM) may 
also b e important if a significant fraction of the ISM is in a cold 
phase jMcKee|[T989l) . A lower velocity (~ 1-2 km s -1 ; lYan et al.l 
12004b achieved in CNM is compensated by a large density (~ 30 
cm -3 ), so that the shattering time-scale in CNM is similar to that in 
WNM (HY09) for graphite. Because of a larger shattering threshold 
velocity of silicate, silicate grains may not be shattered in CNM. 

3.2 Shattering in WIM 

In this Letter, we treat shattering in WIM as an additional effect, 
since larger fraction of the volume is occupied by neutral medium. 
We adopt the results at 100 Myr in Section [3~T1 for the initial con- 
dition. In Fig. [2] we show the grain size distribution after 1 and 3 
Myr of shattering in WIM. By shattering in WIM, not only large 
grains but also relatively small grains with a ~ 0.01-0.1 [im 



are destroyed efficiently. Thus, shattering in WIM may be an ef- 
fective mechanism to enhance the abundance of small grains at 
a < 0.01 \im. This is further discussed in Section |4~T1 in terms 
of the PAH abundance. 



4 DISCUSSION 

4.1 Comparison with the Galactic grain size distributions 

This paper provides an experimental investigation as to whether the 
mean abundance of small grains derived from the mean extinction 
by MRN can be achieved by shattering. Observationally, the grain 
size distribution is derived from the extinction curve along a line of 
sight, which could cross not only diffuse ISM and but also clouds. 
The typical time-scale of the phas e exchange bet ween diffuse ISM 
and clouds is less than 100 Myr dlkeuchil ll988). shorter than the 
grain lifetime. Therefore, the grains processed in WNM are mixed 
among the ISM phases. If this is true, we can directly compare 
the results obtained in this paper with the grain size distribution 
derived observationally from the extinction curve. If we consider 
the existence of the ISM phases other than WNM, the shattering 
time-scale should be modified as T s jmt/ /wnm, where /wnm is the 
fraction of time spent in WNM. We assume that /wnm is not much 
smaller than 1 . 

As shown in the previous section, shattering in WNM redis- 
tributes large grains into small grains in the grain lifetime so that the 
abundance of small grains becomes consistent with the MRN grain 
size distribution. Therefore, shattering in turbulent ISM is a strong 
candidate for the dominant production mechanism of small grains. 
The short shattering time-scale (r s hat ~ 100 Myr) indicates that 
the loss of large grains by shattering is likely to be balanced with 
the dust growt h in dense clo uds, not by the injection of large grains 
from stars (see lDraineT2 009 for the time-scales of these processes). 
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Figure 2. Same as Fig. [T] but with additional shattering in WIM. The shattering durations in WIM are 1 and 3 Myr for the dot-dashed and solid lines, 
respectively. The dashed line shows the initial (i.e. processed by shattering in WNM for 100 Myr; same as the dot-dashed lines in Fig.[TJ, and the dotted line 
presents the MRN size distribution. 



Shattering in WIM destroys not only large grains but also 
relatively small grains very efficiently. Thus, some imprints of 
WIM on the grain size distribution may be seen for the Galac- 
tic grain size distribution. In Fig. [2] we indeed observe the en- 
hancement of small grains at a < 100 A for graphite. The en- 
hancement of very small grains is relatively minor for silicate, 
since silicate is harder to shatter than graphite. It is interesting to 
point out that the enhancement of very small grain abundance is 
more prono unced for carbonaceous speci es th an silicate species ac - 
cording tolWeingartner & Drainej J200lh and lli & Draind d200lh . 
iLi & Draind d200ll) have shown that an excess in the PAH abun- 
dance at a < 100 A relative to the extension of the MRN size 
distribution is required to reproduce the PAH emission strength 
in the Milky Way. Such an excess also keeps the co nsistency 
with the extinction curve 1 Weingartner & Drainel l20o"lh . Strictly 
speaki ng, small PAHs shou l d be t reated as molecules as done by 
iMicelotta. Jones. & Tielens] d2010h . Nevertheless, we expect that 
PAHs once produced are not much affected by shattering in turbu- 
lence, since small grains (also PAHs) are coupled with small scale 
(i.e. low velocity) turbulence. 



4.2 Coexistence of large and small grains 

Our hypothesis investigated in this paper, that is, that shattering 
governs the production of very small grains, naturally explains the 
ubiquitous coexistence of large and small grains. Indeed, in the 
Galactic environment, the spatial distribution of very small grains 
radiating in the MIR and that of large grains emitting in the far in- 
frared (FIR) shows a tight correlation ( Shi bai. Okumura. & Onakal 
1999). Also the PAH emission is always associated with the dust 
continuum dOnaka et alii 19991) . Moreover, the FIR colours of ex- 
tragalactic objec ts are consistent with the Galactic FIR colours 
jHibi et aLl l2006). which implies that the grain size distribution is 



not very different. The fact that observational MI R to FIR SEDs can 
be described by a small numbers of parameters jDale et aT]|200lfc 
iNagata et al]|2002h also implies that there is a ubiquitous mecha- 
nism for determining the grain size distribution. Finally, shattering 
occurs predominantly in diffuse neutral and ionized medium, which 
is consistent with the d iffuse nature of the 70 [im excess coming 
from very small grains ( iBernard et alj2008h . 

We have shown that even if there are only large grains, the 
supply of small grains are efficient enough. In other words, shat- 
tering by interstellar turbulence always provides a path for small 
grain production. After an enough number of small grains are pro- 
duced, it would be more reasonable to consider that the grains sup- 
pled from stars or dense clouds are well mixed with the MRN-like 
size distribution. The shattering under an MRN size distribution 
is considered in HY09, who also obtain a similar time-scale for 
shattering. The similarity in the time-scale between our case and 
HY09's is explained as follows: In WNM, only large grains ac- 
quire velocities large enough for shattering. Thus, a collision with 
a large grain is the only way for a grain to be shattered. Let us 
consider a 'projectile' with size a colliding with a large grain with 
size ai argc . The contribution from projectiles with sizes between a 
and a + da to shattering of the large grain is roughly proportional 
to 7raf arge Khat?i(a)Vda, where V is the shattered volume. If the 
projectile is large, V is large (in particular, if a ~ 0.1 txm, V is 
equal to the entire volume of the large grain). Thus, larger projec- 
tiles have larger impacts on the large grain. Indeed, if we consider 
shattering among large (a > 0.08 p.m) grains only, the production 
rate of grains with a ~ 0.01 |J.m becomes only 30 per cent less (i.e. 
30 per cent of small grain production is attributed to the collisions 
between small and large grains). Since n(a) is similar between our 
initial distribution and the MRN at large sizes, we obtain a similar 
shattering time-scale to HY09. 
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4.3 Metallicity dependence 

Since the shattering time-scale is inversely proportional to the dust- 
to-gas ratio (equation [3}, it is expected that the abundance of very 
small grains is low in dust-poor (meta l -poor) galaxies. However, as 
pointed out by iHirashita & Ichika wa (2009), detection of infrared 
dust emission from low-metallicity objects could be biased to dense 
systems. In such dense environments, shattering is enhanced. Con- 
sequent ly, the MIR radiation from very small grains may become 
strong <Gallianoetal]|2005l : lEngelbracht et alj|2008h . Thus, dust 
emission from diffuse (low-density) low-metallicity system can be 
a strong test for our theory: if the MIR continuum is suppressed 
in such a system, our hypothesis that shattering is the production 
source of very small grains is supported. 

The production of small grains by shattering also indicates 
that the materials of small grains reflect those of original large 
grains. The formation mechanism of large grains should depend 
on metallicity: The production of large grains by accretion of 
heavy elements in molecular clouds is inefficient in low metallicity. 
Thus, stellar production of grains should be the do minant source of 
large grains in low-metallicity environ ments (e.g. iMatsuura et al.l 
l2008MGalliano. Dwek. & Chanialll2008l) . Considering that molecu- 
lar clouds are the site of various chemical reactions, large grains 
formed by stars and those grown in molecular clouds could be 
chemically different. The metallicity dependenc e of MIR emission 
properties, especia lly that of PAH emission (e.g. lSmithetalj2007l : 
iDraine et alj|2007r) . may reflect the different properties of origi- 
nal large grains in terms of metallicity. Also if there is a spatial 
variation of production sources of large grains, spatial variation of 
very small grains an d P AHs may also be pr esent as observed by 
iBernard et all J2008h and lParadis etail b009h . 



5 CONCLUSION 

We have theoretically investigated interstellar shattering of large 
grains (a ~ 0.1 \im) as a production source of small grains. We 
have shown that by shattering in WNM the abundance of small 
grains reaches the level consistent with the MRN size distribution 
in ~ 10 s yr (i.e. in the grain lifetime). Shattering in WIM addition- 
ally destroys grains with a ~ 0.01 |J.m and redistribute them into 
smaller sizes. This effect is more pronounced for the carbonaceous 
grains, and can explain the enhancement of very small carbona- 
ceous grains (or PAHs) wit h a < 100 A as indicated observation- 
ally by lLi& Draine] feOOlri . Since turbulence is ubiquitous in the 
ISM, our theory naturally explains the strong MIR-FIR correlation 
observed generally for galaxies. 



ACKNOWLEDGMENTS 

We thank T. Nozawa and T. Kozasa for helpful discussions, and H. 
Yan for providing the data of grain velocities in turbulence. 



REFERENCES 

Bernard, J.-P, et al. 2008, AJ, 136, 919 
Bianchi, S., & Schneider, R. 2007, MNRAS, 378, 973 
Borkowski, K. J., & Dwek, E. 1995, ApJ, 454, 254 
Dale, D. A., Helou, G., Contursi, A., Silbermann, N. A., & Kol- 
hatkar, S. 2001, ApJ, 549, 215 



Draine, B. T. 2009, in Henning Th., Grim E., Steinacker J., eds., 
ASP Conf. Ser., Vol. 414, Cosmic Dust — Near and Far, ASP, 
San Francisco, 453 

Draine, B. T, et al. 2007, ApJ, 663, 866 

Dwek, E. 1998, ApJ, 501, 643 

Engelbracht, C. W., Rieke, G. H, Gordon, K. D., Smith, J.-D. T, 

Werner, M. W., Moustakas, J., Willmer, C. N. A., & Vanzi, L. 

2008, ApJ, 678, 804 
Gail, H.-P, Keller, R., & Sedlmayr, E. 1984, A&A, 133, 320 
Galliano, F, Dwek, E., & Chanial, P. 2008, ApJ, 672, 214 
Galliano, F, Madden, S. C, Jones, A. P., Wilson, C. D., & 

Bernard, J.-P. 2005, A&A, 434, 867 
Gauger, A., Balega, Y. Y, Irrgang, P., Osterbart, R., & Weigelt, G. 

1999, A&A, 346, 505 
Groenewegen, M. A. T. 1997, A&A, 317, 503 
Hibi, Y, Shibai, H., Kawada, M., Ootsubo, T, & Hirashita, H. 

2006, PASJ, 58, 509 
Hirashita, H., Nozawa, T, Yan, H., & Kozasa, T. 2010, MNRAS, 

in press darXiv: 1 001.2606) 
Hirashita, H, & Ichikawa, T. T. 2009, MNRAS, 396, 500 
Hirashita, H., & Yan, H. 2009, MNRAS, 394, 1061 (HY09) 
Hofmann, K.-H., Balega, Y, Blocker, T, & Weigelt, G. 2001, 

A&A, 379, 529 
Ikeuchi, S. 1988, Fundam. Cosm. Phys., 12, 255 
Jones, A. P., Tielens, A. G. G. M., & Hollenbach, D. J. 1996, ApJ, 

469, 740 

Jones, A. P., Tielens, A. G. G. M., Hollenbach, D. J., & McKee, 

C. F. 1994, ApJ, 433, 797 
Kozasa, T, Hasegawa, H., & Nomoto, K. 1989, ApJ, 344, 325 
Li, A., & Draine, B. T. 2001, ApJ, 554, 778 
Mathis, J. S„ Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 

(MRN) 

Matsuura, M., et al. 2008, in Kwok S., Sandford S., eds., IAU 

Symp. 251, Organic Matter in Space, Cambridge Univ. Press, 

Cambridge, p. 197 
McKee, C. F. 1989, in Allamandola L. J., Tielens, A. G. G. M., 

eds., IAU Symp. 135, Interstellar dust, Kluwer, Dordrecht, p. 43 1 
Micelotta, E. R., Jones, A. P., Tielens, A. G. G. M. 2010, A&A, 

510, 36 

Nagata, H, Shibai, H, Takeuchi, T. T, & Onaka, T. 2002, PASJ, 
54, 695 

Nozawa, T, Kozasa, T, Habe, A., Dwek, E., Umeda, H, Tomi- 
naga, N., Maeda, K, & Nomoto, K. 2007, ApJ, 666, 955 

O'Donnell, E., & Mathis, J. S. 1997, ApJ, 479, 806 

Onaka, T, Mizutani, M., Tomono, D., Shibai, H, Nakagawa, T, 
& Doi, Y. 1999, in Cox P., Kessler M. F. eds., ESA-SP 427, The 
Universe as Seen by ISO, ESA, Paris, p. 731 

Paradis, D., et al. 2009, AJ, 138, 196 

Shibai, H., Okumura, K, & Onaka, T. 1999, in Nakamoto T. ed., 
Star Formation 1999, Nobeyama Radio Observatory, Nobeyama, 
p. 67 

Smith, J. D. T, et al. 2007, ApJ, 656, 770 

Weingartner, J. C, & Draine, B. T. 2001, ApJ, 548, 296 

Yan, H., & Lazarian, A. 2003, ApJ, 592, L33 

Yan, H., Lazarian, A., & Draine, B. T. 2004, ApJ, 616, 895 

This paper has been typeset from a Tr-X/ fflpX file prepared by the 
author. 



